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A BLADE -EIEMENT ANALYSIS FOR LIFTING ROTORS THAT IS 
APPLICABLE FOR LARCffi INFLOW AND BLADE ANGIES AND 
ANY REASONABLE BLADE GEOMETRY 
By Walter Castles, Jr., and Noah C. New 


SUMMARY 


Simple approx ima te solutions are derived for the relationships 
between the rotor thrust and flight-path velocity components and the 
rotor blade angle, torque, and in-plane forces. These approximate solu- 
tions, based upon the assumption of a triangular distribution of blade 
circulation and a parabolic variation of blade-element profile drag with 
lift, are sufficiently accurate for preliminary calculations and the 
determination of the equilibrium angle of attack and lateral tilt of the 
tip-path plane. 

A set of more exact blade-element equations is then derived giving 
the relations between the thrust and flight-path velocity components and 
the equilibrium blade angles, torque, and in-plane forces and moments. 
Neither the blade-element nor the approximate solutions are dependent 
upon the usual approximations that the inflow angle and blade angle of 
the blade elements are small angles. Thus the present equations should 
be useful for convertaplane as well as helicopter calculations. 

It appears that nonlinear blade twist may be desirable for a con- 
vertaplane rotor in order to obtain useful propeller efficiencies. 
Therefore, the blade-element equations have been arranged so that any 
reasonable distribution of blade twist may be used. Also, the equations 
were set up in terms of an arbitrary blade-chord distribution since it 
was found that the use of the actual blade-chord distribution and the 
elimination of the usual assumption that the blade airfoil extended 
inboard to the axis of rotation largely eliminated the necessity for the 
usual reverse-flow corrections. Tables of the necessary factors are 
included for blades having a linear taper, linear twist, and an airfoil 
contour from r = 0.15R to r = R and for blades having a linear taper, 
helical twist, and airfoil contours extending from r = 0.20R to r = R 
(where r is the radius of the blade element and R, the radius of the 
blade tip). 
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The present analysis is based upon the following assumptions: 

(1) The blade-element lift coefficient may be assumed to be pro- 
portional to the sine of the blade-element angle of attack, and the 
blade-element profile-drag coefficient may be represented by the first 
three terms of a Fourier series in the blade-element angle of attack. 

This implies the neglect of blade stall effects in the equations for the 
blade angles. The effect of tip stall is taken into account in the 
equations for the rotor torque. 

(2) The blade axes may be assumed to be, and to remain, straight 
lines . 

(3) The lateral and longitudinal variations of the normal component 
of the induced velocity at the tip-path plane may be assumed to be linear. 

(4) The effects of compressibility on the tip sections of the 
advancing blade may be neglected. 

(5) All radial velocity components and the tangential components 
of the induced velocity may be neglected. 

(6) Blade tip effects may be neglected. 

A comparison of the results given by the present equations with the 
full-scale helicopter test data of NACA TN 1266 shows good agreement for- 
the helicopter flight range covered in .that report. 


INTRODUCTION 


This project, which was conducted at the Georgia Institute of 
Technology Engineering Experiment Station under the sponsorship and with 
the financial assistance of the National Advisory Committee for Aeronautics 
was undertaken in order to develop a blade-element analysis for lifting 
rotors that would be useful for convertaplane calculations. This neces- 
sitated the elimination of the usual approximations that the blade-element 
inflow angle 0 and the blade angle B are small angles and required a 
reasonably exact treatment of the blade geometry. 

It was found that a practical approach to the problem of eliminating 
the small-angle approximations for the lift forces could be obtained by 
writing the lift coefficient of the blade element as 


C Z ~ a sin = a (sin 0 V cos 0 y + cos 0 y sin 0 V ) 
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and, consequently, the thrust component of force dL cos on a blade 
element as 

dL cos 0 V = £ Pa (u cos 0 V ) |sin 0 v (u cos 0 V ) + cos 0 v (u sin 0 v )j c dr 

Similarly, the tangential component of the lift on a blade element may be 
expressed as 

dL sin 0 V = § P a ( u sln 0 v )Jj5in 0 v( u cos 0 V ) + cos e v( u sln $v)\ c dr 


It was also found that the -small- angle approximations could be largely 
eliminated for the profile-drag terms by expressing the blade-element 
profile-drag coefficient c^ as 


PdQ = 6 o + - e l + ^2 cos cLj, 


The exact blade geometry has been retained throughout by expressing 
the blade-chord and blade-twist distribution in the form of the following 
constants : 




dx 


1 r ^ 

a nc=^R c C0S e t xn_1 

'J v-. 


dx 


and 


"ns 


= — i— r c sin 04-x 11- '*' 

rtR J v t 


dx 
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where 6- j. is the blade twist in the angle of zero lift between the 
reference station and nondimens ional radius x. Values of these con- 
stants are given in tables 1 to 3 for blades having linear taper, linear 
twist, and xp = 0.15 and in tables 1, 4, and 5 for blades having linear 

taper, helical twist, and xp = 0.20. 

The present system of equations has been get up with respect to tip- 
path-plane coordinates or coordinates based on the virtual axis of rota- 
tion (fig. l) rather than the usual coordinate system based on the plane 
of zero, feathering in order to obtain shorter expressions for the in-plane 
rotor forces and moments. The use of coordinates alined with the virtual 
axis of rotation also facilitates the treatment of some accelerated flight 
problems . 

Certain refinements in the induced-velocity theory, as given in 
reference 1, have been incorporated with some minor changes in the 
present equations along with the necessary terms for an arbitrary angular 
velocity of roll and pitch of the tip-path plane. 

Standard NACA nomenclature has been used where possible, with the 
subscript v for virtual axis of rotation appended to the usual symbols 
which, in this paper, have a similar meaning but different numerical 
values. 


NOTATION 


a slope of lift curve for blade element at 0.75R (per radian) 

8q rotor coning angle 

a*Q coning angle for zero blade-root bending moment 

ap coefficient of sine component of blade cyclic-pitch angle 

measured with respect to tip-path plane where 

9 V = Aq + - ap sin V + "bp cos i(r 

also coefficient of cosine term of Fourier series for 
blade flapping angle 3 measured with respect to plane 
of zero feathering where 

3 = Sq - ap cos V - bp sin ijr - ^2 cos - ^>2 s ^ n 2ijf - . . . 
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a 2 


A 0 


b 

*1 


b 2 


c 

c o 

c <*o 

c l 

^mx 


Cmy 


coefficient of second-harmonic cosine term in Fourier series 
for blade flapping angle 

mean blade pitch angle at reference station, positive above 
tip-path plane 

number of blades in rotor 

coefficient of cosine component of blade cyclic-pitch angle 
measured with respect to the tip-path plane; also coeffi- 
cient of sine term of Fourier series for blade flapping 
angle measured with respect to plane of zero feathering 

coefficient of second-harmonic sine term in Fourier series 
for blade flapping angle 

blade chord at radius r 

extended blade -root chord at r = 0 (for linear taper) 

section profile-drag coefficient 

section lift coefficient 

rotor rolling-moment coefficient measured about X-axis 

(l^?) 

rotor pitching-moment coefficient measured about Y-axis 



C Q 

AC Qs 


Q 


rotor ‘torque coefficient 
increment to Cq from tip stall on retreating blade 


rotor thrust coefficient ( — 2 


rotor X-force coefficient 


J P* 
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■'xy 


rotor-blade tangential-force coefficient, positive in 

Fxy 


direction of rotation 


p*a 2 B. hl 


rotor Y-force coefficient 


| p^ 2 R 4 j 


L xy 


r y 

F Z 

S 

*1 

X 



C 2 rotor-blade thrust-force coefficient 

I>p fuselage and wing drag 

D 0 blade profile drag 

Eq mean blade drag angle, positive in direction of rotation and 

measured between blade axis and line through rotor axis 
of rotation and drag hinge (i.e., blade drag angle £ is 
£ = Eo + El cos ijf + Fi sin ijr + . . . ) 

Ep coefficient of cosine term in expression for blade drag 

angle 

Fp coefficient of sine term in expression for blade drag angle 

F x component of rotor resultant force acting along X-axis 

tangential component of the resultant air force on blade, 

positive in direction of rotation 

component of rotor resultant force acting along Y-axis 

Z component of resultant air force on blade 

acceleration due to gravity- 

mass moment of inertia of blade about flapping hinge 
nc - 0 nc sin A 0 + 0 ns GOS 


I ns = °ns' sin A 0 " 0 nc cos *0 



MCA TN 2656 


7 


I v mass moment of inertia of rotor about virtual axis of 

rotation 

If mass moment of inertia of blade about drag hinge 

k a blade-root spring constant (blade-root bending moment in 

foot-pounds divided by angular deflection in radians of 
three -quarter-radius point from Sq) 

Lp fuselage and wing lift 


M x 

% 

Q 

r 

r 

r P 

R 

t = 

T 

U 

v 

V 

Vi 

w 


rotor rolling moment 

rotor pitching moment 

rotor torque, negative in direction of rotation 
of blade element c dr 
of blade center of gravity 
of inboard blade airfoil element 
of flapping hinge 
of blade tip 

(for linearly tapered blades) 

rotor thrust, component of rotor resultant force along 
Z-axis 

component of resultant velocity at blade element that is 
normal to blade axis 

mean normal component of induced velocity at tip-path 
plane (positive down and to rear) 

velocity along flight path 

Z component of induced velocity at radius r and azimuth 
angle i|r (positive in plus Z-direction) 

slope of longitudinal variation of nondimen sionad. induced 
velocity • 


radius 

radius 

radius 

radius 

radius 

Tip chord 
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gross weight plus down component of any acceleration force 
acting on aircraft 

nondimensional blade radius (r/R) 

nondimens ional radius outboard of which retreating blade is 
stalled 

nondimensional radius of inboard blade airfoil element 

slope of lateral variation of nondimensional induced velocity 

angle of attack of fuselage measured between flight-path 
velocity vector and longitudinal fuselage axis 

blade-element angle of attack measured from line of zero 
lift 

angle of attack of tip-path plane measured in the XZ-plane 
between flight-path velocity vector and tip-path plane, 
positive below tip-path plane 

blade flapping angle at azimuth angle ijr (for tip-path plane, 
P v = a 0 - a 2 cos 2i]r - b 2 sin 2i | r - . . 
for plane of zero feathering, 

= a 0 " a l cos ^ ~ ^i s * n t - a 2 cos - b 2 sin 2 - . . .) 

circulation of blade element at radius r and azimuth 
angle ijr 

constants in expression for T where r = (Tq + sin i|r)x 

value of c^ at cj = 0 

constant in power equation for cdo 
(i . e . , c^ = 6 q + e c 1 2^ 

constants for first three terms of Fourier series expressing 
relation between c^ and a^. 

(i.e., c^ = 6 o + e l s i n a r + 6 2 cos “r 

or c&q = sin a r + e 2 cos 

blade drag angle at azimuth angle -<{r, positive in direction 
of rotation 
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0]_ twist in zero-lift chord line between axis of rotation and 

blade tip for blades with linear twist, positive for 
increased angle at tip (i.e., 0-fc = 0]x) 

0-^ twist in rotor blade angle of zero lift between reference 

station and radius r, positive for larger angle outboard 

dip design helix angle at tip of blade for blades with helical 

twist 

9 y pitch angle of blade element at radius r and azimuth 

angle ^ measured between zero-lift chord line and tip- 
path plane, positive above tip-path plane 
(A q + % “ a l sin + b-j_ cos ljr) 


e , 


angular displacement of tip-path plane about X-axis from 
horizontal 


e y 

V 

Pv 

P 

ff n 


angular displacement of tip-path plane about Y-axis from 
horizontal 

inflow velocity ratio at center of tip-path plane 
fV sin cty - 


(- 


ftR 




in-plane velocity ratio at tip-path plane 
density of air 


(V cos ay \ 

V Hr ) 




dx constants which express blade-chord distribution 

/ . nl 


i.e. , ^*1 — 


itR , 


c dx 


X 1 


CTo = — / cx dx, etc.) 

>* J x x ) 
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°hc - 


CT ns 


i / c cos e-tx 11 " 1 dx 
Jx 1 

— f c sin O-tx 11 -! dx 
jtR 


constants which express blade-chord and 
blade -twist distribution 


0c 

0V 

“x 

ft 


angle between flight path and horizontal, positive below 
horizontal 

inflow angle at blade element measured in plane perpen- 
dicular to blade axis and between tip-path plane and 
relative wind, positive below tip-path plane 

azimuth angle of blade axis measured about Z-axis from X-axis 
(This iangle is very nearly but not identically equal to the 
equivalent angle in plane of zero feathering. ) 

ratio of angular velocity of roll of tip-path plane about 
X-axis to ft 

ratio of angular velocity of pitch of tip-path plane about 
Y-axis to ft 

mean angular velocity of rotor blade axis about Z-axis 
All angles are in radian measure. 




AHAIYSIS 

Value of Normal Component of Induced Velocity at Radius r 

and Azimuth Angle 


It Is shown in reference 1 that for a lightly loaded single rotor 
composed of' a large number of blades b each having a circulation given 
by the expression 


r = r 0 + r l sin ^ 


( 1 ) 
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the mean value of the normal component of the induced velocity is 



( 2 ) 


Equation (2) was derived on the assumption that the wake extended 
to infinity and had the form of a straight elliptic cylinder. Thus, for 
those flight conditions where a "vortex ring" type flow exists, equa- 
tion (2) is not applicable and the value of v must, at present, be 

obtained from experiment. The term ^1 - ^ p v ^ in the denominator of 

equation (2) arises from the lateral dissymmetry in the blade circulation 
that is required for rolling-moment equilibrium, and this term is the 
only correction which the elementary theory makes in Glauert's original 
hypothesis that v = T/2pAV', where V* is the resultant velocity at 
the center of the rotor. 


If the distribution of the normal component of the induced velocity 
Vi over the tip-path plane is denoted by a power series in the non- 
dimensional radius x and a Fourier series in the azimuth angle such 
that for the first-order terms 


— = + wx cos V + yx sin (3) 

GR OR 

it can be shown from the results of reference 1 that 



y ~ 2p v ^ (5) 

For level flight and p v > 0.15 the expression for y may be simplified 
to 


yxCip 


( 6 ) 
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' It may be noted that for a pair of equally loaded, coaxial, counter- 
rotating rotors, the values of w and y are 


v 


k 

3 



and 



(7) 


y = 0 


(8) 


Approximate Values of Rotor Blade Angles, Torque, and 


X-Force and Y-Force Coefficients 


It is convenient for preliminary calculations and checking and nec- 
essary, in the general -case, for the determination of the angle of attack 
and lateral tilt of the tip-path plane to have simple expressions of use- 
ful accuracy for the rotor torque, X force, and Y force that are 
independent of the rotor blade angles. One such set of equations which 
take into account all the principal variables including the primary 
effects of the reverse-flow region may be obtained from a consideration 
of the distribution of the blade circulation. It may be noted before 
proceeding that the use of any radial blade -circulation distribution 
other than the uniform value assumed in the derivation of the induced- 
velocity equations will underestimate the induced torque. Thus, it is 
theoretically incorrect to calculate the Induced torque from blade- 
element equations. However, for the extreme case of a triangular dis- 
tribution of circulation along the radius and n v = 0-5, the error in 


the induced torque is only 3^- percent and thus is probably within the 
errors of present equations for the induced velocity. 


For the present purposes a triangular distribution of blade circu- 
lation along the radius and a sinusoidal variation with azimuth angle 
is sufficiently accurate and will be used. Then 


r = (fg + r-]_ sin tyjx 


(9) 
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Writing 


U cos = &R(x + M-v sin ''t f ) 


( 10 ) 


and 


U sin 0 V = flR 


Xy + yx sin ijr + (wx - aQ(i v ) cos ijr 


( 11 ) 


it follows for thrust and rolling-moment equilibrium that 


3*f2R 2 C T / \ 

r = — (l - - sin (12) 

b(l - Hv 2 ) V 3 ' 

The reference blade angle Aq at r = 0.75R corresponding to the 
average value of the circulation and inflow angle at this station and 
a weighted chord may be obtained from the substitution 

2ira 0;75R = ^°^0.75R = 2 ( f / cU )o.75r 


or 


Aq = sin 


-1 


irba ^(l - Hy 2 ) 


1 + ~ Hv 2 



+ 2560 

729 




(13) 


where 



cx 11- -'- dx 


(l4) 


For linearly tapered blades the values of cr n may be obtained by 
interpolation from table 1. 



14 


NACA TN 2656 


The values of a^_ and b]_ 

circulation and inflow angle at 
for aj and = 0 and i}r = jt 


obtained from the differences in blade 

r = 0.75R for t(r = it/2 and i|r = 3 rt / 2 
for b^ are 


a l 


y - ^ ^-vPv 8C^n v 

(* - | “vj 3*to 3 (l - f Hv 2 + f Kv 1 *) 


(15) 


and 


bl = _w + - ao^v 


( 16 ) 


where the value of sq for blade root moment equilibrium on a blade 
with the flapping hinge at the axis of rotation is approximately 

3pjcR5c T 6 - - - Pv 2 ) 

ao = (17) 

^Il(l - Pv2) 


Similarly 


a 2 


8aoP v 2 

9 - 8p v 2 


(18) 


and 


b 2 = 0 . ( 19 ) 


The value of the blade-element profile-drag coefficient may be 
represented with sufficient accuracy for the present purposes by two 
terms in a power series in the blade-element lift coefficient c 1 such 
that 


Cdo = B 0 + €C Z 2 


( 20 ) 
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where for conventional airfoils 6q sj 0.0080 and € rj 0.0080. Ma kin g 
the necessary substitutions and integrations the value of the rotor 
torque coefficient is 


Cq = 


c tv~ _ ^v/ 
- dv 2 


+ |' b5 0 


cr I4. + — O^v 2 + dv 2 )^ 



1 + f dv 2 Wl^ 

+ ACqs 

(1 - dv 2 ) 2 


from equation ( 89"b) (where applicable) 


( 21 ) 


Similarly, neglecting blade-shank drag which is assumed to be included 
in the helicopter parasite drag, the values of the X- and Y -force 
coefficients are 


and 


| p*ft 2 R U 

C T (2X^ v - y) _ | ‘^(Isf) , , 

= + b& 0 dv<72 ~ y (22) 

1 - ^ a-nv 2 )2. 



(23) 


The above equations based upon a triangular distribution of blade 
circulation along the radius and a sinusoidal variation with azimuth 
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angle are sufficiently accurate for preliminary calculations, checking, 
and the determination of the angle of attack and lateral tilt of the e 
tip-path plane provided there are no large areas of the rotor outside 
the reverse-flow region that have blade elements operating in the stalled 
or negative thrust range. This implies a reasonable blade twist for the 
flight conditions . 

Table 6 shows a comparison of the values of the parameters 
calculated from the above circulation equations with the flight test 
results of reference 2. 


Determination of Angle of Attack and Lateral Tilt 
of Tip-Path Plane 

Given the values of the flight-path velocity V, climb angle 0 C , 

gross weight and vertical component of the inertia force W, fuselage 
and wing drag, lift, moment characteristics, and position of center of 
gravity, the fuselage angle of attack and thus the fuselage and wing 
lift Lp and drag Dp can be obtained for the trim condition by 
setting the summation of moments, acting on the fuselage and wing and 
taken about the rotor hub, equal to zero. Since the lateral tilt of 
the tip-path plane has a negligible effect, it follows from the geometry 
of the above forces, as shown in figure 2, that 


Dp cos 0 C - Lp sin 0 C + F x cos 0 

tan 0y = 

W - Lp cos 0 C - Dp sin 0c + F x sin 0y 


( 24 ) 


is a good approximation for unaccelerated flight. It may be noted that 
Dp should include an allowance for rotor-hub and blade- shank drag. In 

general, the terms involving F x will have only a small effect on the 
value of 0 y and a sufficiently exact solution can be obtained on the 
second iteration. Thus, as a first approximation. 
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Dp cos 0 C - Lp sin 0 C 
tan Qy — 

W Lp cos ,0c - Dp sin 0 C 


(25) 


Ov = 0c + e y 


(26) 



(27) 


^v 


V COS Oy 

DR 


(28) 




V sin Oy 
DR 


DR 


(29) 


The values of v/DR may he obtained from equation (2) or by double 
interpolation from table 7 which includes the experimental values for 
vertical descent from reference 3 and estimates of the values for the 
inclined flight "vortex ring" states. The values of w, y, and F x 

can then be determined from equations (4), (5), and (22), and from these 
second approximations to the values of Qy, cty , . and n v can be made 

from equations (2k), (2 6), and (28). If necessary, a new value of Cip 
may then be obtained from the equation 


Ct£ — 


W - Lp cos 0 C - Dp sin 0 C + F x sin Qy 


P«D 


v; 


cos 0, 


(30) 


and thus the more exact value of ' Xy from equation (29). 


For helicopter calculations the first approximation for Gp is 
sufficiently accurate, and if is small (i.e., (j. v < 0.15) the effect 

of F x on ay may be neglected for level flight. ‘ 
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The tail-rotor thrust Trp required for a helicopter with a single 
main rotor is 


T t = Q /l 


(31) 


where l is the perpendicular distance' between axis of main and tail 
rotors and the value of Cq may he obtained from equation (2l) • The 

lateral tilt 8 X of the tip-path plane for a single-rotor aircraft in 
unaccelerated flight is thus 


2 C y + C Q f 


'X ~ 


Cfp 


where C v is given hy equation (23). 


(32) 


Application of Two-Dimensional Airfoil Theory and Data 
to Rotor-Blade-Element Calculations 
Two-dimensional thin-airfoil theory demonstrates that 


cj = a sin a (33) 

For a two-dimensional cascade of airfoils, equation (33) is modified 
hy a multiplying function of the solidity, chord spacing, and hlade angles 
that is very nearly unity for average lifting- rotor configurations as 
shown in reference 4. Thus, within the approximation that the radial 
components of flow may he neglected, equation (33) should he applicable 
for blade-element rotor theory over the unstalled range of blade-element 
angles of attack. Beyond the stall, equation (33) is somewhat less in 
etror than the usual relation cj = aa as can he seen from figure 3 

which is a plot of the above expressions and the experimental values of 
Cj against a for an NACA 0015 airfoil. The use of equation (33) > 

rather than the usual approximation that c j = aa, allows the thrust and 

tangential components of lift on a hlade element to he exactly expressed, 
within the approximations involved in neglecting radial components of 
the flow, in terms of the easily integrated in-plane and normal components 
of the velocity at the hlade element U cos 0 V and U sin 0 V . Thus the 

usual approximation that the inflow angle 0 y is a small angle may he 
eliminated. This may be demonstrated as follows: 
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Omitting the negligible component of the profile drag, the thrust dT 
on a blade element c dr is 


dT = i pU 2 cc-£ cos dr 


( 3 ^) 


or since 


= a 


sin a^, = a (sin 0 y cos 0 y + cos 0 y sin 0 y ^ 


(35) 


dT = ^ pac(u cos 0 y ) jjsin 0 y (u cos 0 y ) + cos 0 y (u sin 0 y jj dr ( 36 ) 

The tangential component of the lift on a blade element may be similarly 
expressed as 


dL sin f v - ^ pac(u sin 0 y ) sin 0 v (u cos 0 y ) + cos 0 y (u sin 0 y ) dr (37) 


The value of the slope of the lift curve a of the blade -element air- 
foil in the above relations may be taken as the value corresponding to 
the Reynolds number, Mach number, and surface roughness existing at the 
three-quarter-radius point of the rotor blades under consideration. For 
the usual tip speeds, in the 500-foot-per-second range, the Prandtl- 
Glauert Mach number correction , 


^Jl - M 2 


(38) 


where 

a' low Mach number lift-curve slope from two-dimensional 

wind-tunnel tests 

M free-stream Mach number at three-quarter blade 

radius 

may be used to correct the lift-curve slope from low Mach number data. 


The values of c^ obtained from two-dimensional wind-tunnel tests 
at appropriate Reynolds numbers and model surface roughness should be 
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directly applicable to rotor-blade- element calculations in the unstalled 
range of angles of attack below the Mach numbers and angles of attack 
for drag divergence, since the effect of subsonic Mach number on profile 
drag is negligible as shown in reference 5* However, it should be noted 
that the profile-drag coefficient is only constant with change of sub- 
sonic Mach number if it is taken as a function of the lift coefficient. 
If the airfoil section data are plotted against section angle of attack 
and the Prandtl-Glauert correction is applied to the lift- curve slope, 
this is equivalent to multiplying the section-angle-of-attack scale by 
o 1/2 

(l - M 2 ;' . Consequently, the section-angle-of-attack scale on the 

profile-drag curve must be multiplied by (l - to retain the 

same relation between c^ and C2- 


In view of the errors in the magnitude and distribution of the blade 
circulation that arise from the necessary neglect of blade deflections, 
and so forth, it is probably not Justifiable to take into account sec- 
ondary effects of the profile drag. Thus, expressing the relation between 
the profile-drag coefficient and the blade-element angle of attack by the 
first three terms of a Fourier series gives 


c do = e 0 + e l sin V'r + e 2 cOS 


(39) 


The constants in the above equation may be evaluated from the two- 
dimensional wind-tunnel data for the blade airfoil at, say, a = 0°, 5°, 
and 10°. The advantages of equation (39) over the usual expression 


c do = 5 o + 5 l a r + 5 2 a r 2 


are: The last two terms of equation (39) can be exactly expressed in 

the known velocity components U cos 0 V and U sin 0 V ; the resulting 

expressions for the forces and moments on the blade are considerably 
simplified by the absence of the squared term in c^; and it is an 

equally accurate approximation to the experimental values of c^ as 

may be seen from figure 4. However, in using equation (39) it may be 
noted that the calculated value of c^ is the small difference between 

large quantities and thus the values of 6 q, e and e 2 should be 
determined to four places in order to obtain the value of c^ to the 
customary accuracy. For the more severe convertaplane flight conditions 
where the inflow velocity is large (|Xy| > O.io) a certain error arises 
in the treatment of the ?q terms, and it is necessary to fall back on 
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the two-term approximation for c^, = 6 i sin ct-p + e 2 cos 

where and e 2 are evaluated from, the experimental data at, say, 

a = 2° and a = 7°. This additional’ approximation is permissible for 
these flight conditions, since the relative effects of the profile drag 
become less important as the inflow velocities and rotor blade angles 
increase. For example, in propeller calculations the single-point 
approximation c^ = €C2 is usually used. 

It follows from the geometry and equations (35) and (39) that the 
tangential component of the profile drag on a blade element -may be 
expressed as 


dD 0 cos 


0 V = \ P c ( u cOS 0v)|^ 6 O u + 

e l(( u cos $v) sin e v + ( u sin 0v) cos e v ] + 

e 2 [(u cos 0 V ) cos 0 V - (u sin 0 V ) sin 0^J“ dr (4o) 


The thrust 



Thrust of a Blade at Azimuth Angle ifr 
F z of a blade at azimuth angle iff is 

c^U cos cos 0 V ) sin 0 V + (u-sin 0 V ) 


cos 0, r dr 


v] 


(41) 


where r^ is the radius of the inboard blade airfoil element. In the 
general case it follows from the geometry that 

!J cos = £ffi(x + p v sin iff) (42) 


and 


U sin 0 V = 2R jly + (wx + nyx - a Q n v ) cos ijr + (y - cd x )x sin iff + 


2b2X cos 2 ^ - 2a2X sin 2i^J 


(43) 
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where 

oijj ratio of angular velocity of roll of tip-path plane about 

X-axis to ft 

ay ratio of angular velocity of pitch of tip-path plane about 

J Y-axis to ft 

Neglect ing the higher, harmonics of the cyclic pitch that may arise from 
control- system linkages, the pitch angle 0 V of a blade element at 
radius r and aximuth angle ijr, measured with respect to the tip-path 
plane , is 


6 V = Aq + 9^ - aj sin i|r + b^ cos i{r 


m 


where 

Aq mean blade pitch angle at reference station 

twist in rotor blade angle of zero lift between reference 
station and radius r 

minus the coefficient of sine component of blade cyclic- 
pitch angle measured with respect to tip-path plane 

b^ coefficient of cosine component of cyclic-pitch angle 

measured with respect to tip-path plane 

In the general case (i.e. , for the convertaplane ) Aq and 6^ 

may not be small angles'. However, it appears that the magnitude of the 
cyclic-pitch angle will always be limited by tip stall on the retreating 
blade to the range where it is a good approximation that 


sin 


^-a-^ sin i}r + bj 


cos 



-a-j_ sin + b^ cos i|r 


m 


and 


cos 


^-a^ sin i)r + bj cos 
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It follows from equation ( 44) , upon expanding the functions sin 0 V 
and cos 0 V , that 

sin 0 V = [sin A q + (cos Ag)^-a 1 .sin i|r + bj_ cos i{r[J cos 9 ^ + 

[cos Aq - (sin A^|(-a-j_ sin ij( + cos \{r[j sin 0^ (A 7 ) 


cos 0 V = [os Aq - (sin sin ijr + bj_ cos T|r[| cos 0^ 


[in Aq + (cos A^^a^ sin + b^ cos i|rjj sin 6 ^ 


m 


Substituting the values of U cos U sin 0 V , sin 0 V , and cos 0 V 

from equations (A 2 ), (A3), (A7), and (¥)) in equation (Al), defining 


J nc 


itR 


f! ^ 


cos 8 ^ dx 


(A 9) 


J ns 


= ^ T cx I ^” 1 sin 0 t dx 

Jx i 


(50) 


I nc = CT nc s in A 0 + u ns cos Aq 


(51) 


Ins - c ns sin Aq - cr nc cos Aq 


( 52 ) 


multiplying out the terms, and reducing the functions of tJt to 


harmonic form give for the thrust coefficient C~ = 


Z \ prt^R^ 

blade at an azimuth angle i)r the expression of equation (53): 


of one 



I 3 c 

i + | a i(y - < D x) - 


1 . / \ 
g °UV + wyf 




-£ a i(y-®*)- 

coa 2? 

5 tl(v + my) 

sin 3^ 



* 2 c 

*lc 

i(a 0 + * 2 )^ 1 ^’.' ' 

5 eibgi T 

£( a lV + Hv)l^v 

( 2 - + a lK + 

l a i(y - ®*W - 
$ M w + ®y)^v 

£ SflVv 2 

a 0 ^ 2 * 1 v “ b lV + 

\ a l( w + ®y)v^ - 
\ bl(y - ®t)mt 

\ 

~t a O a l' 1 y 2 ‘ 

-(1 a 0 + a 2 ) a lhr + 

hVr 

(| a O tl l + " 

2 V> 


-$ a i(y -.®> v - 

t M w + ®yK 

fr °oVv 2 

a i( v + “vK + 
l bi(y - 

4 v " T 


ro 

P- 


Equation (53) 


J 3e 

hs 

lie 

^( e i a 2 + t i 1 ^)( v + “y) + 

1 /n.V^ n Vs _ \ / lr _ ,*i\ \ 

- ^d^L) \* ~ 

*lUv ■ Vr " 

J(y - <%)n v - 

^( a l a 2 + b 1 b 2 )e 0 n v 


»i - (y - “it) 

(nybg - Bjby)^ + 
b2M v 

(j[ a y j v - K)»v 

-by - (w + tBy) 

( a o + + 

(Oydg + bybg)^ 

- 1 b lMy £ 

2a 2 + a l®2(y “ ^fx) - 
a 2 b l( w + “V) 

( a 0*2 - OVy " 
K v + “VK 

| a 0Mv 2 

- 2 b 2 - ejb 2 (y - m*) + 

V2(* + “*) 

\ 

-( a l + o 0 b l b 2)l J v + 

|(y - %K 




a l ^ 



^ bm., 2 
4 ■*■■ * 
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Equation (53) is written in tabular form where the coefficients in the 
boxes must be multiplied by row and column heads. Values of oh c and 
^ may be obtained by interpolation from tables 2 and 3 for linearly 
tapered and twisted blades , where 


c = c Q (l + tx) from x = 0.15 to x = 1 (54) 


6 t = djx 


(55) 



(56) 


and 

c 0 extended blade-root chord at r = 0 

t _ Tip chord 
c o 

twist in angle of blade, zero lift between axis of rotation 
and tip 

In order to use the tabulated values of cr nc and cr ns for blades with 

linear twist and taper, it is necessary to take the reference blade 
pitch angle Aq at the extended blade -root chord c Q at r = x = 0. 

The use of the lower limit x^ = 0.15 in the computations for the 

blades having linear taper and twist corresponds to present practice 
and largely eliminates the necessity of making any reverse-flow correc- 
tion to the blade thrust. The reverse-flow effects' are discussed in the 
following section. 

Additional tables (tables 4 and 5) give the values of. cr nc and a ns 
for blades having linear taper from x-[_ = 0.20 to x = 1 and helical 
twist where 

^tan 0rji^ 


0^ = tan' 


(57) 
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and is the design helix angle at x = 1. In this case, the refer- 

ence station for A q is taken at the blade tip. The tables for helical 
twist are included for convertaplane usage since helical twist would 
appear to he desirable for a reasonable propeller efficiency. An inner 
limit of Xj_ = 0.20 was used for the computation of the values of 
cr nc and cr^ for this case of helical twist in order to minimize the 
severe root stall likely to occur under some convertaplane flight con- 
ditions. It might be pointed out that helical twist would also appear 
to afford an increase in helicopter-rotor performance over that obtain- 
able with linear twist. 


Reverse -Flow Considerations 

For normal helicopter and convertaplane flight conditions 'where 
there is a downflow through the rotor and 0 V is negative over the 
reverse-flow region, the maximum value of p v is limited for conven- 
tional rotors to relatively low values of the order of 0.30 by tip stall 
on the retreating blades. Under these conditions the portion of the 
retreating blade extending inboard from the outer edge of the reverse- 
flow region at x = -p v sin ijr, where the in-plane component of velocity 

is zero, to x = x^, where the blade airfoil section ends, has a neg- 
ligible thrust loading because the in-plane components of velocity are 
very small. The present equations take into account the fact that the 
blade airfoil does not exist inboard of x = x^, for which region the 
the in-plane components of velocity are larger, within the reverse-flow 
circle, and previous equations erred in assuming the blade airfoil to 
exist. 


For those flight conditions where there is an upflow through the 
rotor and the tip-stall limitations on p v are relaxed, the present 

equations give the proper direction to the blade -element thrust for 
those blade elements within the reverse-flow region and inside the radius 
where j v x 20 v . 

Thus, for all practical purposes, it is not necessary to use reverse- 
flow corrections when applying the present equations to conventional rotors. 

For unconventional rotors operating with net downflow at large values 
of p. v it would appear from strip analysis to be desirable and even nec- 
essary to minimize the forces in the reverse-flow region by using a suf- 
ficiently large design value of x^, for example, > ~ ^v max ' 

this case less error is introduced by taking into account the inboard 
blade airfoil limit and omitting the usual reverse-flow correction than 
vice versa. 
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Mean Rotor Thrust 

Omitting the coefficients of the second- harmonic flapping angle 
which have a negligible effect on the mean rotor thrust, the value of 
the mean rotor thrust coefficient obtained from the first row of equa- 
tion (53) by averaging the value of C z /a over the interval from 
i|r = 0 to if = 2it is given by 

Ur = C 1 + \ a l( y - ®x) - \ b l( v + °v)] T 3c + I a 0 b lM2c + 


|( a lV + Hv)Mlc - \ - a l^ + |(y - ®xVv J 2s (58) 


Mean Rotor Air Rolling Moment 

The value of the mean rotor air-rolling-moment coefficient about 
the X-axis 


C 


mx 


Mx 

Tj' pTttt 2 R 5 


is found, upon integration, to be obtained by multiplying the second row 
of equation (53) by ^ b and introducing the moment arm by writing the 
subscripts of I nc and I ns to one higher order. Thus, 


2C 


mx 


ab 


[^v + a lV + | a i(x - - 1 b l( w + -y)^3c + 

i a 0 b 1 p v 2 I 2c + (a x - y + a*)!]*, + 2 s 


(59) 


Mean Rotor Air Pitching Moment 


Similarly, the mean rotor air-pitching-moment coefficient 


C my 


7j P«2 2 r5 
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obtained from the third row of equation (53) is 


^ = [W " l a l( w + + | b l(y " %)^v) I 3c + 

| a 0 a lt i v 2 I 2c + ^ + w + oy)!^ - a 0 ^ v I 3s + ^ b^ 2 !^ 


(60) 


Mean Blade-Root Air Moment 

The coefficient of the blade-root air moment Mq is merely 
the first row of equation (53) with the I factors to one higher sub- 
script. Thus, for 


, Mq 
J m ° ' | 


^° = |l + | ai (y - | b x (w + ay)] Ii^ + i a 0 b l^v I 3c + 

i(a!A. v + P- v )^vl2 c - \\ ~ a l^v + ^v] I 3s 


Equilibrium Values of Mean Rotor Pitching Moment 
and Rolling Moment \ 

If an exter nal moment M]_ is applied to a single rotor with three 
or more blades about a diameter, axis 1, the differential equations of 
motion about axis 1 at V = ti and- axis 2 at t = tp + 90° can be shown 
by the use of Euler' s equations to be 


1 dab 

1 — i + ftCJ Up + 

2 dt 2 


^l^L = % 
I v Iv 


(62) 
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and 


dCDo koCDp 

— - - no>! + ■ - = 0 

dt x I 


(63) 


where cd^ and <xq are the angular velocities of the tip-path plane 
about axes 1 and 2, respectively, k^co^ and ls>b>> are the damping 
moments, and 1^. is the mass moment of inertia of the^ rotor about the 

virtual axis of rotation. The general solution of equations (62) 
and ( 63 ) is a pair of equations of the form 


CD. 


lor 2 



2 

t + B cos 




k l +k 2 

Iv 


t 


(6k) 


In the actual case, damping of the nutation appears to be very rapid 
for an articulated rotor. Also, for pilot-controlled motion, k 2 ~ 0. 

For example, for a constant control moment Mj, k 2 = 0, and k^ = 22I V , 
which is then the value of k^ for critical damping, it follows that 


or 


CD]_ = 


2M 1 t 


e -2flt 


(65) 



(66) 


It can be seen from equations ( 65 ) and (66) that the transients decay 
very rapidly and their effects can be neglected in most problems. There- 
fore, to a good approximation for a single rotor 


“x = + “rf 


(67) 


% 


= -IvO 2 




( 68 ) 
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where and My-p are any moments transmitted about the X- and Y-axes 

from the fuselage to the rotor. For steady straight and level flight 


.tOjj = ay = 0 


■ (69) 


For steady hanked turns the value of Cp can he taken proportional to 
sec 0x • Also 


-g sin 9v tan By 

<% 2 — - (70) 

vn 

and 




g sin 9 X tan 9 2 

vn 


(71) 


where 9 X is the equilibrium lateral-tilt angle of the tip-path plane 
(approximately equal to equilibrium angle of bank, positive for turns 
in direction of rotor rotation) . 

For any curvature of the flight path, the nondimens ional compon- 
ents qx and ay of the spatial angular velocity of the aircraft may 
be calculated and, consequently, the approximate equilibrium values of 
M x and My can be obtained from equations ( 67 ) and (68). 


Approximate Solution for Equilibrium Values of Mean Reference 
Blade Angle Aq, Lateral and Longitudinal Components of 
Cyclic Pitch ap and bp, and Coning Angle sq 

An approximate solution of the set of four nonlinear, transcendental 
equations ( 58 ), (59).> ( 60 ), and (6l) for the four unknowns Aq, ap, sq, 
and bp that is sufficiently accurate for most steady-flight helicopter 
work and useful as a first trial for steady-flight convertaplane calcula- 
tions may be obtained as follows: Setting the small terms and <x>x, ay, 

and Cmx equal to zero and cos Aq = 1 in equations (58) and (59) and 
eliminating ap gives 


sin Aq 


(1? " CT 3s “ M2c)( ff 4c + J Pv 2ff 2c - V3s) + ^v 2(T 2c( 2ct 3s + V2 c) 
( a 3c + § Pv 2 < 7 1 c - M2s)( ct 4c + 4 ^v S<J 2c ~ \°3b) ~ %hr 2 °2c a 3c 
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Then, from equation (59) for ay = ay = C^. = 0 


a l 


^v^c ^4s _ ^rhylgs 
(v + £ ST l v) i 3c + I 4s + f ^v 2 i 2 s 


(73) 


Let a 0 he the design coning angle for the general case of semirigid 
hlades (i.e., coning angle for zero blade-root bending moment). Let 
k Q Q be the spring constant of the blade for angular deflections of the 

three- quarter-radius point from Hq. Then setting the summation of 
moments about the blade root equal to zero and solving for aQ, the 
coning angle at the three-quarter- radius point, 

\ PJtan^jlllc + §(*l*v + ^v)M2c ~ a l^v) I 3s] + ^0% " ^g 


1 + 


(74) 


where 


mass of blade 

radius of blade center of gravity 

mass moment of inertia of blade about flapping hinge 
(or root) 

If the blades have a flapping hinge at the axis of rotation a Q = k S Q = 0. 
If the flapping hinge is located at radius r p from the axis of rotation. 


a n = 0 and kg. 


- (r p /0.75R) 


Then, knowing a Q , it follows from 


equation (60) that for 


c% - ay - C^y - 0 


*1 


3 oM3e 


wl 

T 


4s 


(75) 


\^3c + I 4s + ], ^v -*-2s 
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For those steady, unaccelerated flight conditions where cos Aq z; 1 , the 
above solutions are sufficiently accurate and may be used to calculate 
the blade loadings and rotor torque, X force, and Y force. 


"Exact" Solution for Aq, aq, and b]_ for Accelerated Flight 
Conditions and Those Flight Conditions where cos Aq ^ 1 


A reasonably rapid and sufficiently accurate solution of the "exact" 
equilibrium equations given by the first three rows of equation ( 53 ) can 
be obtained by using an approximate value for the coning angle sq such 
as that given by equation (IT) or (7^) • 


Then for the approximate value of Aq given by equation (72) and, 
for example, two other values several degrees successively smaller, the 
"exact" corresponding values of aq and bq can be determined by 


rewriting the equilibrium equations for the rotor pitching and rolling 
moments in the form 

2 C„ 


Aai + Bb^ — P 


my 


and 


ab 

2C„ 


Ca n + Db, = R + — SSL 
11 ab 


(76) 


where 


A - |(w + co jr )^ T I 3c - | aQp v 2 I 2c 


(77) 


® “\r^3c ]j.(y °)x) ix v-*-3c " ^^s “ • 1 j, Pv^2s (78) 

f 

\ 

C - ^v^3c + a )x) 1I v^3c + T ks + ^ Pv^2s (79) 

D = - I(w + a V )p v I 3c + i aQp v 2 I 2c (80) 


P = (w + oy)l 4s - a 0 p v I 3s 

R = -2p v I 3c + (y - oj 1^ + Vv I 2s 


(81) 


(82) 
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Then 


and 
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Having computed the values of a-^ and bj for each of the assumed 
values of A Q , the corresponding values of Gj, may he found from the 
equation for the thrust equilibrium where 


2&P 

ah 


1 2 i 

+ 2 I l c ~ ' >c v I 2s 2^ y ai x)^v I 2s + 


[g( y u ^c)- I -3c + 2 ^v^v-^lc + ^v^-2s| a l + 

[| a 0M2c - |( w + a V) I 3c] b l (85) 

Then plotting the values of 2CT/ab, ap, and bp against the trial values 
of Aq } the "exact" value of Aq, and thus ap and bp, may he obtained 
from the plot at the design or desired value of Crp. 
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In-Plane Component of Force F^y on a Blade 
at Azimuth Angle i|f 

The in- plane component of force in the direction of rotation F-jjy 
on a blade at azimuth angle is from equations (37) (4o) 


= | pa J c(u sin 0 V ) J^sin 8 v (u cos 0 V ) + cos 8 v (u sin 0 v jJ dr 


l p / r c(o cos dv){« 0 U + e^Jsin 0 v (u cos 0 ir ) + 


0v) 


cos 0 v (lJ sin 0 V ^J + e 2 |cos 0 v (tF cos 0 V ) - sin 0 v (u sin 0 V ) 


■dr 

( 86 ) 


where 


c do = 6 0 + 6 1 sin a r + e 2 cos “r 


Then, by (l) substituting the previously evaluated expressions for 
U cos 0 V , U sin 0 V , sin 0 V , and cos 0 V given by equations (42), (43), 
(47), and (48); (2) neglecting the effects of second-harmonic flapping; 
and (3) writing (coU)(u cos 0 V ) as e 0 (u cos 0 v ) 2 (u/u cos 0 v ) and 
expanding 



by the binomial theorem and dropping third and higher terms, the expres- 
sion for the constant and first -harmonic terms becomes 


Cxy " “ (^)a - - (fiCzy) ei - (&xy )< 2 (87) 



X 3 C 


Bln + 

y - ^ ai(y - a ^) 2 - 

\ b i(y - “*) ( v + »r) + 

i a i( v + °v ) 2 

COS f 

V + CDy - ^ b 1 (y - CD *) 2 + 

| »i(y - %)(v + <*y) - 
^ bi^v + ay ) 2 


I2c 

He 

\ ^ 1/v » V, ^ 

,w y T " '-‘-xj^v * 

»i(y - “x)^ - 

b l( w + «Jr)*v 


*• a o b i(y - ^xW - 

| * “y)^ 

( Q 1^V + t»v)\r + 

j a o 2s i> l ir 2 

_a 0* x v “ 

1 aofti(y - ®je)hv + 
1 a 0 b l( w + “V) 1 ^ 

- b lV 2 - 

| a 0 S^ v 2 


































| e O b lJ i T §(? 3 >v + Mt)m v 

+ a l*v + 

l «l(y - “r>y - I “oVv 2 

£ bi(v + my)^ 

^ “l( w + “j) Hy ' 

l bl(y - ®x)n v - 


-t *o“i ^ 2 
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Rotor Torque 

The effects of tip stall at the higher values of p v and C^/cr^ 
on Cq are large and may he approximately evaluated for high-speed 
flight where Xv is negative, as follows: The retreating blade will 

he stalled outboard of the nondimens ional radius (for X v negative) 


X v + p v tan 
Xs = 

y + tan^SSi - Aq - a x - A0 

where A 0-^ is the aerodynamic blade twist between the reference station 
and the tip. Assuming a jump of 0.08 in the value Of c^ at the stall 

and that the rotor area within which blade stall exists is a' segment of 
minimum radius x s and symmetric about i|r = 3^/2, the increment ACq s 
to Cq due to tip stall is approximately 



hoax 


- Aq - a i - A0 t J 


( 8Qn 3 


ACq s « 



- p v ) 2 (l - 


x s ) \[T^ 



( 89 b) 


(if x s < gv or x s > 1 equation ( 89 b) is not applicable and 
ZC Q =0.) Then 


2Cq , 

— S s -(Constant terms of Cyy with subscripts n on a n , I nc , and 
b J 


ns 


increased to n + l) + 


2ac 


Qs 


(90a) 


For steady-state calculations equation (90a) may be reduced to 
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2 Cq _ ^ ~ M + ^Qs + 

b b(l - p v 2 ) 13 

e Q J^l + i w 2 )ff4 - i a o wn v 0 3 + |(^ 2 + ^v 2 )^ + 

eijl^c +| a 0 t ll i v I 3c '+ |(n v + a 1 X v )n v I 2 c + ( a l^v - ^v^s] + 

6 2 RaiUv “ x v) I 3c - - \ ^l^) 1 ^ ~ \ a o' b l^v I 3 s “ 

2 ( a l^v + t x v)( a v 1 2sj 

Rotor X Force 

The value of the rotor X-force coefficient C x is 

= -^Sine terms of Cxy^ ( 9 l) 

However, the greater part of C x arising from the lateral variation in 
blade circulation is a small difference between large quantities which 
are principally functions of a^_ and Aq • It follows that this part of 
C x is more accurately obtained from the circulation equations than from 
the blade-element equations. - Thus for steady- state solutions 


2C X 

b 


2Ct(2A. v u v - y) 


+ e 0 (2ji v + y^v) a 2 + 


b(l - Hv 2 ) 

^ij^l^v + 2p v )l2c + ( a l ~ y) I 3s + ^v “ ^vj^vlis| + 
e 2^ a l “ y) I 3c + a l^v “ ^v^M-v^lc " (2dv + a l^v)^2s] (92) 


Rotor Y Force 

The value of the rotor Y-force coefficient is 


2 C V 

= Cosine terms of Cj^y 


( 93 ) 
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As in the expression for C x ' the above value of Cy given by the blade- 

element equation is a small difference between large quantities and the 
result for steady-state flight is more accurately obtained from the 
circulation expression 


C T (w - | aQUy) 
1 - Pv 2 


(94) 


Second-Harmonic Flapping 

Again letting ka^ be the spring constant relating the blade-root 

bending moment in foot-pounds to the angular deflection in radians of the 
three-quarter-radius point of the blades from the unstressed position, 
it follows that the magnitude of the cosine component of the second 
harmonic of the blade flapping angle is 


a2 


J + KL 
1 - KM 


Similarly the magnitude of the sine component is 


(95) 


b 2 


^ L + JM 
~ 1 - KM 


(96) 


where 

J 


(terms not involving b 2 in the cos row of thrust equa- 
tion (53) with the I factors changed to one higher sub- 

\ l -k 

script) x ( — 

\3Iifl 2 - kaoj 


K 


(coefficients of b 2 in the cos row of thrust equa- 
tion (53) with the I factors changed to one higher 

\ / ^ pitn 2 R5 a 

subscript) v ‘ d 


x 


(31 1^ - k. 


*0; 


(terms not involving a 2 in the sin 2f row of thrust 
equation (53) with the I factors . changed to one higher 

\ / J P*l2R5 a 

subscript) x I — 

\3Ii2 2 - kaoy 
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itO 


M 


(coefficients of a 2 in the sin 2ijr row of thrust 

equation (53) with the I factors changed to one higher 
v / i pjtfl 2 R5 a 
subscript^ X | OT 


31^^ - kao / 


For steady-state flight conditions where ay = ay = 0 the expres- 
sions for the factors 3, K, L , and M may he simplified to 


i pita 

3Ixn 2 - ka 0 


[j-|( a l^v + ^)M2 c " ( a l " 


(97) 


K: 


i pitaf2 2 R5 


SI]? 2 - k 


a 0 


(- 2i 4b) 


m 


i pitafl 2 R5 
3Il« 2 - k a 


0 


I2 \ b l X v^ i v I 2c " (^1 + 2 V ) |a v I 3 s + 2 a 0^v 2 l 2s 


(99) 


M 


^ pjtaH 2 R5 

3iiH 2 - k Bg 



( 100 ) 


and 1-^ is the mass moment of inertia of the blade about the flapping 
hinge . 

It may be noted that kg 0 = 0 for, blades having a flapping hinge 
at the axis of rotation. If the flapping hinge is located at radius rp 
then 


w n2 




a 0 


1 - 


'A. 


0.75R 
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Amplitude of Constant and First-Harmonic Components of 

Lag Angles in Unaccelerated Flight 

For an articulated rotor having lag hinges normal to the plane of 
rotation and located at a small radius e the equilibrium blade lag 
angle Eq is 


E, 


§ P*r 5 r 2C ~1 - 

» from equation ( 90 ) 

M(.e fl - — — \ L * J 

5 V o.7R/ 


(101) 


where is the mass moment of the blade about the lag hinge. 

Similarly the coefficients of the cosine and sine components of the 
lag angle are 


En 


and 


PrtR^ - 2a 0 b l6 l£ 

M^ e - h 


■g- PrtR^F ^ + Sa^a-i R I^ 

M^e - 1^ 


( 102 ) 


(103) 


where a-^ s and b-^ B are the cos and sin components of the angle 

between the tip-path plane and the hub plane. For unaccelerated flight 
the values of a^ s and b-j_ s are approximately 


*16 * °Sr - “f 


(104) 


-’Is 


®xf " ®x 


(105) 


xf 


where 6. 


is the equilibrium lateral tilt of the fuselage 
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Also 






mass moment of inertia of a blade about lag hinge 

coefficient of cos in equation (88) for with sub- 

scripts of I factors changed from n to n + 1; approxi- 
mate value from circulation equations is 


= 


2c t(| w - ^v) 

*(l - ^v 2 ) 


(106) 


coefficient of 


sin in equation (88) for C^y with sub- 


scripts of I factors changed from n to n '+ 1; approxi- 
mate value from circulation equations is 


f £ = 


2Ct 


— X-yM-vJ 


h(l - Pv2) 


0 •O85Crp 2 p v 0^ 

0 . 008 ( 2 n v + y\v)ao + — — , . 

b 2 Q3 2 (l - p v 2 ) 2 

(107) 


Thrust Unbalance 

Two-bladed rotor .- The second -harmonic variation in Gp for a two- 
bladed rotor is 


ACip 

a 


= Fourth + fifth rows of equation (53) 


(108) 


For (% = o>y = 0 and steady-state conditions, the equation for the 
amplitude may be simplified to 


ACrp 

a 


|[f a 2 I 3s - ( b l + I »)M2 s + 5 V/ 1 J + 

, \2l l/2 

( 2 b 2 I 3 s + a l^ x v I 2sj r 


( 109 ) 
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Three -b laded rotor .- The third-harmonic variation in Op for a 
three -h laded rotor is approximately 


2A&P ' - 

s; Sixth + seventh rows of equation (53) 

3a 

or for tpx = ay = 0 the amplitude is approximately 


(110) 


2ACy 

3a 





( 111 ) 


An Independence -of -Blade -Element Analysis for Hovering, 

Vertical Ascent , and Convertaplane Propeller Condition 

The use of the relation c^ = a sin a permits a considerable 
' simplification of the equations resulting from the assumption of the 
independence of blade elements. As the exact propeller solutions of 
Betz, Goldstein, and Theodorsen are not applicable to a lifting rotor 
at zero or small advance ratios, a simple analysis of the independence 
of blade elements may be viseful. 

From momentum considerations the thrust dT on an annulus of the 
rotor disk 2«r dr is related to the induced velocity Vp at the rotor 
element by the expression 


WTS ■ Vi ( Vi + 7 8111 “*) (112) 

/ 

But 


+ V sin 0 ^ = 11 sin 0 V 


(113) 


Thus 


— (u 

4jtpr dr 


sin 


0v) ( u sin 


V sin 



( 11*0 



hh 
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The thrust of the annulus is also equal to the thrust acting on the 
portions of the blades -within the annulus which is 


dT = ^ pbU^ccj cos 0 V dr (115) 

where 

c-£ = a sin Oj. = a ^sin 0 V cos 0 V + cos 0 V sin 0 v j (ll 6 ) 

Thus 

dT = — pab^U cos 0 V ) jsin 0 v (u cos 0 V ) + cos © v (u sin 0 v )j c dr (117) 


Substituting the above values of dT in equation (ll4) and solving for 
U sin 0 V 


U sin 0 V /Vo aba r \ 

-(f + ^r cose i - 


HR 


fv a abcr r 


_ + _ c ° s e v j + 


1 2 aba T 


q x sin e v 


(118) 


where 


v a = 


V sin a v 

Hr 


(119) 


cr,, = 


itR 


( 120 ) 


Then from equation (117) 


2G 

ab 


T 


-It 


x sin 0 ir + 


sin 0 , 


HR 




cos £? v UpX dx 


( 121 ) 
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U sin fiy 

where the value of — at -x is given by equation (ll8)'. 

2R 

Similarly, from blade-element considerations 


TT - - a / ( ” m4) [i s±n e v + ( u 5 g- 0v ) cos a* 


r 1 c do r /U sin 0 _\ ' P 

/ — : x sin 0 V + -) cos Q„\ajcsr dx 

J Xl sin a r L v \ OR ) Y J r 


( 122 ) 


where the values of 


'do 


are obtained from a plot of 




sin a r ^ sin a r 

a r for the blade airfoil at values of a r given by the relation 


against 


0^ 


= © v + tan' 


-1 


1/ lJ sin 0 V \ 
x\ flR / 


( 123 ) 


If it is necessary to take into account the rotation of the slip- 
stream for large rates of vertical ascent or the propeller condition, 
this may be accomplished to a first approximation by using an effec- 
tive ft, ft e , in every case where 



( 124 ) 


The geometry of the above equations. is exact and they are convenient 
for graphical or numerical integration on account of the repetition of 
factors . 

Neglecting the induced radial and tangential velocity components, 
the optimum blade-angle distribution for minimum induced power and a 
given blade-chord distribution and nondimens ional axial flight-path 

, U sin (ZL 

velocity v a may be obtained by setting — - equal to the constant 

value giving 



h6 
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sin 9 V - 


^v(^v “ v a) 




1 

[V 2 - * 2 

)\f- 

[V - v a)5] 

( 

\ ~ v a] 

2x2 J 


where 


n"h rr 



The optimum chord distribution for a given desired constant value of 
along the "blade and the same restrictions is 


8^v(^v - "V"a) 

°r = . 

"bc^^v 2 + x 2 

For this optimum chord distribution,, the optimum distribution of 6 y 
reduces to 


sin By = 



1 4 - \ 

, , a 2 V 2 - c l 2 (K 2 + x2 ) 

\ 


L \ 

c i x J 


For calculations where the flight-path velocity and equilibrium 
value of Crji are known or can be estimated, the following procedure 
may be followed: 

(1) Calculate and plot the radial distribution of u r 

(2) Calculate the effective value of Crp and v a where 



(125) 

(126) 

(127) 
c l 

(128) 

(129) 
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(3) Calculate the approximate value of Aq 
which for these flight conditions reduces to 


from equation (72) 


/2CT, 


sin A 


0 


\T5 Ef - °3s - V*2c, 
( ff 3c " V CT 2s) 


(4) Calculate and plot the radial distrihution of 0 V = Aq + 0-^ 
for the value of Aq obtained under item (3) and two lower values at 
increments of several degrees 

(5) Calculate and plot the radial distribution of U sin 0 v /fl e R 
for the above distribution of 0 y from equation (ll8) using (i = fl e 
throughout 

(6) Calculate and plot the radial distribution of the integrand 
of equation (121) for the three values of* Aq and graphically or 

numerically integrate for the values of 2(>p e y/ab corresponding to the 
three values of Aq 

(7) Obtain the correct, value of A 0 from a plot of C^ against Aq 

(8) Calculate and plot the radial distribution of the integrand of 
equation (122) for the three values of Aq and graphically or numerically 

Integrate for the values of SCQ^b corresponding to the three values 

of A 0 


(9) Obtain -the equilibrium value of C^ at the equilibrium value 
of Aq from a plot of CQg against Aq 


(10) Calculate the equilibrium value of 



Comparison of Experimental and Calculated Values of Parameters 

Table 6 shows a comparison of the experimental data of reference 2 
for those runs where Gj ~ 0.005^-5 with the values calculated by the 

approximate blade-element equations of this report. The blade-element 
lift-curve slope was assumed by the authors to have been a = 6.5 from 
the experimental results of reference 6. The values of gq., 
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and ^2 were evaluated for the points = 0.0090, 0.0105, and 0.0170 

at ci = 0, 0.5, and 1.0, respectively, from figure 19 of reference 6. 

The exact solutions for the various parameters differ from the 
tabulated approximate solutions by a negligible amount for these heli- 
copter flight conditions. 

A consideration of the results presented in table 6 would 
indicate that much of the remaining discrepancy between experimental 
and calc ula ted blade angles and torque coefficients may be due to the 
neglect, in the present calculations, of the effects of the rotor induced 
velocity on the lift and drag of the fuselage. 

It may be noted that the longitudinal component of the angle 
tan - l/— 1 between the rotor resultant force and the thrust component 

V 2C T ) 

normal to the tip-path plane is very small for all these helicopter 
flight conditions and that the direction of the resultant is inclined 
forward for those flight conditions where there is a net downflow through 
the rotor. The inclinations of the tip-path plane to the horizontal 0 X 
and 0y are also small angles and, consequently, for many unaccelerated- 
flight helicopter calculations the rotor resultant force can be assumed 
to be perpendicular to the tip-path plane and the thrust equal to the 
gross weight without introducing serious errors. 


CONCLUDING DISCUSSION 

' l 

Simple relations for the rotor blade angles and the values of Cq, 
C x , and Cy, derived upon the assumption of a triangular distribution 
of blade-element circulation along the radius and a sinusoidal variation 
with azimuth angle 'in conjunction with a linear variation of profile 
drag with lift, would appear to be useful for helicopter and converta- 
plane performance estimation and the determination of the equilibrium 
angle of attack and lateral tilt of the tip-path plane. 

The blade-element equations, based upon the relation .that 
cj = a sin a r = a(sin 0 V cos + cos 0v sin 0 V ) , and the cr nc and 
a ns functions of the blade-chord and blade-twist distribution afford 
a reasonably exact and concise treatment of the geometry and should be 
useful for convertaplane as well as helicopter calculations. 
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The use of the empirical relation + 6 1 s 4n a r + e 2 cos a r , 

rather than the usual expression that c&q = 5q + 5-j_a r + 62a . r ^, consid- 
erably simplifies the equations for the in-plane forces and moments and 
presents a sufficiently exact solution of the geometry for helicopter 
calculations . 

For convertaplane calculations, the approximation that 

c do = e i sin Oj, + ^2 cos ctj. allows an exact treatment of the geometry 

and should be a sufficiently accurate expression for c^ at the larger 

advance ratios where the effects of the profile drag become of less 
relative importance. 

The larger sources of the remaining errors in the blade-element 
analysis probably have the following order of importance for contemporary 
helicopters : 

(1) The neglect of the effects of blade-element stall implied in 

the relation that cj = a sin ctp 

(2) The neglect of the effects of blade flexibility 

(3) The neglect of the radial variation in the normal component 

of the induced velocity 

- (4) The neglect of the effects of compressibility on the tip 

sections of the advancing blade. 

t 

Georgia Institute of Technology 
Atlanta, Ga., May 15, 1951 
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TABLE 1 

VALUES OF a FOR BLADES WITH LINEAR TAPER 
n 


ri 


^Interpolate for values for given t; 


"0 



t = - 1; c = c ft (l + tx) 

c 0 0 ' J 


t 

0 

b 

1 — 1 
to 

a 2/ ff o 

a 3 / a o 

a b / a 0 

X-J_ = 0.15 

0 

0.8500 

0.14888 

0.3322 

0.2499 

-1 

.3612 

.1566 

.0823 

.0499 

= 0.20 

0 

0.8000 

o.48oo 

0.3307 

0.2496 

-1 

.3200 

.1493 

.0811 

.0497 
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TABLE 2 


VALUES OP 0 nc FOR BLADES WITH LINEAR TAPER, LINEAR 


TWIST, AND = 0.15 


Interpolate for values for given t first and then for values for 

^ c n 

given 0-, ; reference station for A n at x = 0; cr n = — ; 

x u u «R 


c tin 


- 1 ; c = c Q '(l + tx); 0 t = 0-^xj 
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TABLE 3 

VALUES OF cr ns FOR BLADES WITH LINEAR TAPER, LINEAR 
TWIST, AND x^ = 0.15 

Jlnterpolate for values for given t first and then for values for 

Q 

given 0 -, ; reference station for A a at x = 0; cr n = -2-; 

x . u u itR 


t = - 1; c = c Q (l + tx); 0 t = 0-jxJ 
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TABLE 4 

VALUES OF 0 nc FOR BLADES WITH LINEAR TAPER, HELICAL 
TWIST, AND x x = 0.20 


Interpolate for values for given t first and then for values of 

c 

given 0mj reference station for A_ at blade tip; a n = — 


jffi 


t = 


'tip 


- 1; c = c Q (l + tx); 0 t = tan _1 ^i tan 0,^ 



ff 2c/ 0 O 


r 3c/°b 


T 4c/ a o 



0.8000 

.7906 

•7654 

.7305 

.6907 

.6489 

.6065 

.5645 

.5231 


0.3200 

.3144 

.3002 

.28o4 

.2594 

.2385 

.2184 

.1994 

.1815 


0.4800 

.4762 

.4651 

.4500 

.4313 

.4io4 

.3882 

.3651 

.3416 



0.3307 0.0810 0.2496 o.o497 


.1016 

.0934 


.3287 
.3233 
• 3149 
.3042 
.2919 
.2782 
.2635 

.2481 


.0803 

.0782 

.0751 

.0714 

.0674. 

.0632 

.0590 

.0547 


.2484 

.2451 

.2398 

.2328 

.2244 

.2151 

.2046 

.1935 


.0493 

.0483 

.0468 

.0448 

.0427 

.o4o4 

.0379 

.0354 
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TABLE 5 

VALUES OF 0 ns FOR BLADES WITH LINEAR TAPER, HELICAL 


TWIST, AND X! = 0.20 


Interpolate for values for given t first and then for values for 

c n 

given 0 q ij reference station for Aq at blade tip; 0 Q = 


t = 


'tip 


- 1; c = c 0 (l + tx) ; 9 + = tan 


'Xl tan ftp) 




















TABUS 6 

cchpabisoj or expbhimskeal akd cALcauum value of pabaheiebs tot 

TBOSB ELKS OF HETIHEKOS 1 FOR WHICH C* s 0.003*3 


VJ1 

a\ 


[All angles In deg) tip stall on run 5j 


Parameter 

Bun 7 

Ct 0 0.005*3) OR - *43 ft/eoc 
Level flight at *3.7 aph 

Eun * 

Ct - 0.00538| OR - **7 ft/ sac 
Laval flight at 58.6 aph 

Run 2 

Ct = 0.005*5 ) or » 445 ft/ 
Level flight at 71.7 *rjfc 











Brparlaantal 

Calculated 

Calculated 

Bxperiamrtal 

Calculated 

Calculated 

Eijarirantnl 

Calculated 

Calculated 



(a) 

<*> 


(a) 

(b) 


(a) 

(t) 

°f 

-a. z 



-*.3 



-6.9 



°f 

23.* 



23.7 



24.1 



<k 

0 



0 



0 



“T 

-2.06 

-1.99 


-3.83 

-3-75 


-3.02 

-5.75 


*0 

7.11 

7 - 1 * 

7.35 

8.17 

7.63 

7-95 

10.10 

' 72 

9-19 

»i 

2.92 

3.08 

2.86 

’ 4.37 

*.*l 

3.84 

6.08 

6.63 

5.05 

■■o 

5.16 

8.2 8 

7.99 

8.30 

8.25 

8.02 

8.67 

8.39 

S.21 

*1 

3.2* 

2.7* 

2.78 

3.30 

3-01 

3.06 

3-93 

3.39 

3.40 

*2 

0.2* 

O.lS 

0.17 > 

0.35 

0.26 

0.26 

0.46 

O.43 

o. 4 e 

ta 

0.00 

0.00 

-0.0* 

-0.08 

0.00 

-0.09 

-0.11 

0.00 

nO .17 

<*> 

0.000202 

0.00020* 

0.000205 

0.0002** 

0.000230 

0.000218 

0.0003*2 

^.oooss* 

a o. 000 328 

10 

-TT45 

- 7.68 

-TI 73 

- 3.83 

-0.03 

-B.lS 



-ji. 56 

% 

0.4l 

0.32 


0.54 

0.37 


0.67 

0 .66 


■ **i 

-0.21 

-0.08 


-0.18 

-0.11 


-0.27 

-O.lB 


c* 


-0.0000*1 

-0.0000*6 


-0.000036 

-0.000062 


- 0.000092 

-o.ooouo 

% 


-0.000263 



-0.000319 



-0.000370 


T “' l fe) 


-o.ci 



-0.30 



-0.48 


6 X 


-0.08 



-0.09 



0.06 



-2.08 

-1.99 


-3-83 

-3.74 


-3.82 

-5.75 - 




0.013 



0.029 



0.058 


\°5 ^ 3 rd 












“Calculated free circulation aquations. 

^Jalculated from blade-eleneent equation*. 
c 3 SQ It added to fuselage drag area f from rafaranc* 7 to 
ollov for rotor-hub blftdft-sbuak, counter -torque- rotor, «ad engiaa- coo Ling dng. 
^In elud e s co i i eobi oc for tip stall. 




®N*cbmn±cal Inpxt subtracted. 
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TABI£ 6 - Concluded 

COMEAHIBOT OT KFEBIMHITAL ADD CALCUMTED VALUES OF PARAMETERS FOR 
1B0SE RUHS OF HEFIRHKE 1 FOB WHICH Cl «u 0.00949 - Concluded 


Op “ 

Hun 11 

0.00946; OR = 443 ft/sec 

525 ft/* ia cliMb at 51*3 

Szporlsazrtal 

Calculated 

Calculated 


(a) 

(1) 


Hen 19 

Cj - 0.00949? OB - 443 ft/aac 
126b ft/eln auto rotative descant at 37-7 mph 


Eiperlmontal 




Calculated from circulation aquations. 

^Calculated frea blade-eloaent aquation*. 

°3 eq ft added to fuselage drag area f fro® reference 7 to allo w 
for rotor-bub, blade-shank, counter-torque-rotor, and. engins -cooling drag. 
e Mecbanical input' subtracted. 


Calculated 

00 
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TABLE 7 


VALUES OF X± = ^ j- FOB GIVEN VALUES OF \ 



AND ^ 


V sin Oy-fe - 3l y 2 
flR f C51 



*>1.28 

*1.15 

*> 1.02 


0 

-.20 

-. 1(0 

-.60 

-.80 

- 1.00 

- 1.20 

-1.40 

-1.60 

-1.80 

- 2.00 

-2.40 

-2.80 

-3.20 

-3.60 

-4.00 

-5.00 

- 6.00 

- 8.00 

- 10.00 


a Experimental . 
^Estimated . 


0.60 

0.80 

1.00 

*0.580 

0.481 

0.457 

*.68 

.543 

.509. 

*>.82 

.630 

• 574 

*>1.03 

•767 

.659 

*>1.42 

1.000 

.769 

*1.77 

1.220 

.896 

*1.65 

1.25 

-97 6 

b l.52 

1.21 

1.000 

*1.39 

1.15 

.984 

*1.27 

1.07 

•947 

*>l.l6 

1.00 

•897 

*>1.06 
13 .96 

•924 

CVJ 

CO 

.914 

.854 

.786 

.833 

.787 

.731 

.765 

.724 

.680 

.699 

.668 

.630 

.64o 

.615 

.586 

.588 

.569 

.544 

.543 

.526 

.506 

.501 

.488 

.472 

.464 

.454 

.44o 

.432 

.424 

.4ii 

.403 

.395 

.386 

.355 

.350 

• 342 

.316 

• 311 

.306 

.282 

.280 

.276 

.256 

.254 

.251 

.234 

.233 

.230 

.192 

.191 

.189 

.162 

.161 

.160 

.123 

.122 

.122 

.100 

.100 

.099 

1 


1.20 1.4o 1.60 1.80 
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TABLE 7 - Concluded 


VALUES OF Xj_ = -I/----- 3 — - FOR GIVEN VALUES OF ^ = Hyl 2 — 3Mv 
fiR| Crji C_ 


V sin a- r , 2 -■ 3m/ 

AND X z = — \ - Concluded 

Z DR 




Z 

2.40 

2.80 

3.20 

3.60 

4 .oo 

5.00 

6.00 

8.00 

10.00 

2 . 4 o 

0.315 

0.285 

0.261 

0.239 

0.210 

0.184 

0.156 

0.120 

0.097 

2.20 

.329 

.295 

.267 

.245 

.224 

.186 

• .158 

.121 

.098 

2.00 

.344 

• 305 

.275 

.250 

.228 

.188 

.159 

.122 

.099 

1.80 

.357 

• 315 

.282 

.256 

.233 

.191 

.161 

.122 

.099 

1.60 

.370 

.325 

.289 

.260 

.237 

.192 

.162 

.123 

.099 

i. 4 o 

.384 

• 333 

.295 

.265 

.240 

.195 

.163 

.124 

.099 

1.20 

.395 

.341 

.301 

.269 

.243 

.196 

.164 

.124 

.099 

1.00 

. 4 o 4 

.347 

.306 

.272 

.246 

.197 

.165 

.124 

.100 

.80 

.413 

,352 

.309 

.276 

.248 

.198 

.166 

.125 

.100 

.60 

.415 

• 356 

.313 

.277 

.249 

.199 

.166 

.125 

.100 

.40 

.416 

• 357 

.312 

.278 

.250 

.200 

.167 

.125 

.100 

.20 

, 4 l 4 

.357 

.312 

.278 

.250 

.200 

.167 

.125 

.100 

0 

. 4 io 

• 354 

.310 

.278 

.250 

.200 

.167 

.125 

.100 

-.20 

. 4 o 4 

.350 

.309 

.275 

.248 

• 199 

.166 

.125 

.100 

-.40 

.395 

.345 

• 305 

.273 

.247 

.198 

.166 

.125 

.100 

-.60 

.386 

• 339 

.301 

.270 

.245 

.197 

.165 

.125 

.100 

-.80 

•374 

.331 

.296 

.267 

.242 

.196 

.165 

.125 

.100 

-1.00 

.362 

• 323 

.290 

.262 

.239 

.194 

.164 

.124 

.099 

-1.20 

.349 

.314 

.284 

.258 

• 235 

.192 

.163 

.124 

.099 

- 1 . 4 o 

.337 

.305 

.277 

.252 

.231 

.190 

.161 

.123 

.099 

-1.60 

.325 

.296 

.270 

.247 

.227 

.188 

.160 

.122 

.099 

-1.80 

.312 

.286 

.263 

.242 

.223 

.186" 

.158 

.121 

.098 

-2.00 

.300 

.277 

.255 

.236 

.219 

.183 

.157 

.121 

.098 

- 2 . 4 o 

.278 

• 259 

.241 

.224 

.209 

.178 

.153 

.119 

.097 

-2.80 

.258 

.242 

.227 

.213 

.200 

.172 

.149 

.117 

.096 

-3-20 

.239 

.226 

.214 

.202 

.191 

.166 

.145 

.115 

.095 

-3.60 

.222 

.216 

.201 

.191 

.182 

.160 

.l 4 i 

.113 

.094 

- 4.00 

.207 

.198 

- .189 

.181 

.173 

.154 

.137 

.111 

.093 

-5.00 

.175 

.170 

.164 

.159 

.155 

.139 

.127 

.105 

.089 

-6.00 

.152 

.148 

.144 

.l 4 o 

• 137 

.126 

.117 

.094 

.085 

-8.00 

.118 

.120 

.115 

.112 

.ill 

.105 

• 093 

.088 

.078 

-10.00 

.096 

.095 

.094 

.093 

.092 

.089 

.085 

.078 

.070 















